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Abstract In this paper, a fermionic hierarchical model is defined, inspired by the Kondo model, which 
describes a 1-dimensional lattice gas of spin-1/2 electrons interacting with a spin-1/2 impurity. This 
model is proved to be exactly solvable, and is shown to exhibit a Kondo effeet, i.e. that, if the interaction 
between the impurity and the electrons is antiferromagnetic, then the magnetic susceptibility of the 
impurity is finite in the 0-temperature limit, whereas it diverges if the interaction is ferromagnetic. 
Such an effect is therefore inherently non-perturbative. This difficulty is overcome by using the exact 
solvability of the model, which follows both from its fermionic and hierarchical nature. 

Keywords Renormalization group • Non-perturbative renormalization • Kondo effect • Fermionic 
hierarchical model • Quantum field theory 


1 Introduction 

Although at high temperature the resistivity of most metals is an increasing function of the tempera¬ 
ture, experiments carried out since the early century have shown that in metals containing trace 

amounts of magnetic impurities (i.e. copper polluted by iron), the resistivity has a minimum at a small 
but positive temperature, below which the resistivity decreases as the temperature increases. One in¬ 
teresting aspect of such a phenomenon, is its strong non-perturbative nature: it has been measured in 
samples of copper with iron impurities at a concentration as small as 0.0005% [14], which raises the 
question of how such a minute perturbation can produce such an effect. Kondo introduced a toy model 
in 1964, see Eq.(2.1) below, to understand such a phenomenon, and computed electronic scattering 
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amplitudes at third order in the Born approximation scheme [13], and found that the effect may stem 
from an antiferromagnetic coupling between the impurities (called “localized spins” in [13]) and the 
electrons in the metal. The existence of such a coupling had been proposed by Anderson [2], 

Kondo’s theory attracted great attention and its scaling properties and connection to lH Coulomb 
gases were understood [11,3,5] ^ when in a seminal paper, published in 1975 [20], Wilson addressed and 
solved the problem by constructing a sequence of Hamiltonians that adequately represent the system 
on ever increasing length scales. Using ideas from his formulation of the renormalization group, Wilson 
showed, by a combination of numerical and perturbative methods, that only few (three) terms in each 
Hamiltonian, need to be studied in order to account for the Kondo effect (or rather, a related effect 
on the magnetic susceptibility of the impurities, see below). 

The non-perturbative nature of the effect manifests itself in Wilson’s formalism by the presence of 
a non-trivial fixed point in the renormalization group flow, at which the corresponding effective theory 
behaves in a way that is qualitatively different from the non-interacting one. Wilson has studied the 
system around the non-trivial fixed point by perturbative expansions, but the intermediate regime 
(in which perturbation theory breaks down) was studied by numerical methods. In fact, when using 
renormalization group techniques to study systems with non-trivial fixed points, oftentimes one cannot 
treat non-perturbative regimes analytically. The hierarchical Kondo model, which will be discussed 
below, is an exception to this rule: indeed, we will show that the physical properties of the model can 
be obtained by iterating an explicit map, computed analytically, and called the beta function, whereas, 
in the current state of the art, the beta function for the full (non-hierarchical) Kondo model can only 
be computed numerically. 

In this paper, we present a hierarchical version of the Kondo model, whose renormalization group 
flow equations can be written out exactly, with no need for perturbative methods, and show that 
the flow admits a non-trivial fixed point. In this model, the transition from the fixed point can be 
studied by iterating an explicit map, which allows us to compute reliable numerical values for the 
Kondo temperature, that is the temperature at which the Kondo effect emerges, which is related to the 
number of iterations required to reach the non-trivial fixed point from the trivial one. This temperature 
has been found to obey the expected scaling relations, as predicted in [ 20 ]. 

It is worth noting that the Kondo model (or rather a linearized continuum version of it) was shown 
to be exactly solvable by Andrei [ 6 ] at h = 0, as well as at h 7 ^ 0, [7], using Bethe Ansatz, who proved 
the existence of a Kondo effect in that model. The aim of the present work is to show how the Kondo 
effect can be understood as coming from a non-trivial fixed point in a renormalization group analysis 
(in the context of a hierarchical model) rather than a proof of the existence of the Kondo effect, which 
has already been carried out in Ref. [ 6 ,7]. 


2 Kondo model and main results 

Consider a 1-dimensional Fermi gas of spin-1/2 “electrons”, and a spin-1/2 fermionic “impurity” with 
no interactions. It is well known that: 

( 1 ) the magnetic susceptibility of the impurity diverges as /3 = —)■ 00 while 

( 2 ) both the total susceptibility per particle of the electron gas {i.e. the response to a field acting on 
the whole sample) [ 12 ] and the susceptibility to a magnetic field acting on a single lattice site of the 
chain {i.e. the response to a field localized on a site, say at 0 ) are finite at zero temperature (see remark 
( 1 ) in App.G for a discussion of the second claim). 

The question that will be addressed in this work is whether a small coupling of the impurity 
fermion with the electron gas can change this behavior, that is whether the susceptibility of the 
impurity interacting with the electrons diverges or not. To that end we will study a model inspired by 


^ The obstacle to a complete understanding of the model (with Ao < 0) being what would later be called the 
growth of a relevant coupling. 
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the Kondo Hamiltonian which, expressed in second quantized form, is 

Ljl-l /A \ 

Ho= X] ( X] -l\ c-{x)^ 

aeftd} a;=-L/2 ^ ^ 

Hk = Ho + Vq + H^ + V ( 2 . 1 ) 

V)) = -AoX] c+(0)o';Li.a2c”,(0)d+3CT^3_„/"^ 

J = 1 . 2,3 
01,0:2,0:3,04 

Vh = -h u>j Y dt^Yda,, 

3 = 1 , 2,3 (a.a')e{Td}" 

where Ao, h are the interaction and magnetic field strengths and 

(1) c^{x),d^, a =t,i are creation and annihilation operators corresponding respectively to electrons 
and the impurity 

(2) CT^ , j = 1, 2, 3, are the Pauli matrices 

(3) X is on the unit lattice and —L/2, L/2 are identified (periodic boundary) 

(4) /\/(x) = f(x + 1) — 2/(x) + f{x — 1) is the discrete Laplacian. 

(5) u> = (wi, u; 2 ,<^ 3 ) is a norm-1 vector which specifies the direction of the magnetic field. 

(6) the —1 term in Hq is the chemical potential, set to —1 (half-filling) for convenience. 

The model Eq.(2.1) differs from the original Kondo model in which the interaction was 

3 

-^oY ,«2 (0) 

i=i 

where is the j-th Pauli matrix and acts on the spin of the impurity. The two models are closely 
related and equivalent for our purposes (see App.A). The technical advantage of the model Eq.(2.1), 
is that it allows us set up the problem via a functional integral to exploit fully the remark that “since 
the Kondo problem of the magnetic impurity treats only a single-point impurity, the question reduces 
to a sum over paths in only one (“time”) dimension” [5]. The formulation in Eq.(2.1) was introduced 
in [6]. 

The model will be said to exhibit a Kondo effect if, no matter how small the coupling Ag is, as 
long as it is antiferromagnetic {i.e. Ag < 0), the susceptibility remains finite and positive as f3 ^ 00 
and continuous os /i —)■ 0, while it diverges in presence of a ferromagnetic (i.e. Ag > 0) coupling. The 
soluble model in [6] and Wilson’s version of the model in Eq.(2.1) do exhibit the Kondo effect. 

Remark: In the present work, the Kondo effect is defined as an effect on the susceptibility of the 
impurity, and not on the resistivity of the electrons of the chain, which, we recall, was Kondo’s original 
motivation [13]. The reason for this is that the magnetic susceptibility of the impurity is easier to 
compute than the resistivity of the chain, but still exhibits a non-trivial effect, as discussed by Wilson 
[ 20 ]. 

Here the same questions will be studied in a hierarchical model defined below. The interest of this 
model is that various observables can be computed by iterating a map, which is explicitly computed and 
called the “beta function”, involving few (nine) variables, called “running couplings”. The possibility 
of computing the beta function exactly for general fermionic hierarchical models has been noticed and 
used in [10]. 

Remark: The hierarchical Kondo model will not be an approximation of Eq.(2.1). It is a model that 
illustrates a simple mechanism for the control of the growth of relevant operators in a theory exhibiting 
a Kondo effect. 

The reason why the Kondo effect is not easy to understand is that it is an intrinsically non- 
perturbative effect, in that the impurity susceptibility in the interacting model is qualitatively different 
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from its non-interacting counterpart. In the sense of the renormalization group it exhibits several 
“relevant”, “marginal” and “irrelevant” running couplings: this makes any naive perturbative approach 
hopeless because all couplings become large {i.e. at least of 0(1)) at large scale, no matter how small the 
interaction is, as long as Aq < 0, and thus leave the perturbative regime. It is among the simplest cases 
in which asymptotic freedom does not occur. Using the fact that the beta function of the hierarchical 
model can be computed exactly, its non-perturbative regime can easily be investigated. 

In the sections below, we will define the hierarchical Kondo model and show numerical evidence 
for the following claims (in principle, such claims could be proved using computer-assisted methods, 
though, since the numerical results are very clear and stable, it may not be worth the trouble). 

If the interactions between the electron spins and the impurity are antiferromagnetic {i.e. Aq < 0 
in our notations), then 

(1) The existence of a Kondo effect can be proved in spite of the lack of asymptotic freedom and formal 
growth of the effective Hamiltonian away from the trivial fixed point, because the beta function can be 
computed exactly (in particular non-pertubatively). 

(2) In addition, there exists an inverse temperature Pk = called the Kondo inverse temperature, 

such that the Kondo effect manifests itself for j3 > Pk. Asymptotically as Aq —t 0, n/f (Aq) = ci|Ao|“^ -I- 

0 ( 1 ). 

(3) It will appear that perturbation theory can only work to describe properties measurable up to 

a length scale in which n 2 (Ao) depends on the coupling Aq between the impurity and the 

electron chain and, asymptotically as Aq —> 0, n 2 (Ao) = C 2 log |Ao|“^ -f 0(1) for some C 2 > 0; at larger 
scales perturbation theory breaks down and the evolution of the running couplings is controlled by a 
non-trivial fixed point (which can be computed exactly). 

(4) Denoting the magnetic field by /i, if > 0 and ^ 1, the flow of the running couplings tends 
to a trivial fixed point (h-independent but different from 0) which is reached on a scale r{h) which, 
asymptotically as /i —)■ 0, is r{h) = Cr \ogh~^ + 0(1). 

The picture is completely different in the ferromagnetic case, in which the susceptibility diverges at 
zero temperature and the flow of the running couplings is not controlled by the non trivial fixed point. 

Remark: Unlike in the model studied by Wilson [20], the T = 0 nontrivial fixed point is not infinite 
in the hierarchical Kondo model: this shows that the Kondo effect can, in some models, be somewhat 
subtler than a rigid locking of the impurity spin with an electron spin[15]. 

Technically this is one of the few cases in which functional integration for fermionic fields is con¬ 
trolled by a non-trivial fixed point and can be performed rigorously and applied to a concrete problem. 

Remark: (1) It is worth stressing that in a system consisting of two classical spins with coupling Aq the 
susceptibility at 0 field is 4/3(1 -|- hence it vanishes at T = 0 in the antiferromagnetic case 

and diverges in the ferromagnetic and in the free case. Therefore this simple model does not exhibit a 
Kondo effect. 

(2) In the exactly solvable XY model, which can be shown to be equivalent to a spin-less analogue of 
Eq.(2.1), the susceptibility can be shown to diverge in the /3 —>■ oo limit, see App.G, H (at least for 
some boundary conditions). Therefore this model does not exhibit a Kondo effect either. 


3 Functional integration in the Kondo model 


In [20], Wilson studies the Kondo problem using renormalization group techniques in a Hamiltonian 
context. In the present work, our aim is to reproduce, in a simpler model, analogous results using a 
formalism based on functional integrals. 

In this section, we give a rapid review of the functional integral formalism we will use, following [8, 
17]. We will not attempt to reproduce all technical details, since it will merely be used as an inspiration 
for the definition of the hierarchical model in section 4. 

We introduce an extra dimension, called imaginary time, and define new creation and annihilation 
operators: 


c±(x, tf=^ e‘^«c±(x)e-‘^», dt{tf= 


(3.1) 
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(3.2) 


for a G to which we associate anti-commuting Grassmann variables: 

Ca ^ V’a {x,t), (t) I—^ (t). 

Functional integrals are expressed as “Gaussian integrals” over the Grassmann variables:^ 

f Pid^)Pidij) ■ /’n^(^¥^«)Wa) • (3.3) 

a 

P{d(pa) and P{d'ipo,) are Gaussian measures whose covariance (also called propagator) is defined by 


r Tre l^^^c^{x,t)c+{x',t') . 


/ / I if\ j 

g-4,,a{X - X ,t-t) = < 


Tr e-/5ffo 


if t > F 


Tre ^^“c+(x',t')c„ (x,f) . 


(3.4) 


Tr e-/5^o 


if t < F 


9v,a\t t ) — 


[' Tid^{t)d'^{t') a t > t' 
-Trd'^{t')d~{t) if f < F 


dsf 

By a direct computation [8], Eq.(2.7), we find that in the limit L,/3 —>■ oo, if e{k) = (1 — cosk) — 1 = 
— cos fc (assuming the Fermi level is set to 1, i.e. the Fermi momentum to ±^) then 




dkodk 

-iko + e{k) ’ 


dko ) 

27r —iko 


(3.5) 


If /3, L are finite, J in Eq.(3.5) has to be understood as ^ '^ko Z ^fc’ is the “Matsubara 

momentum” kg = ^ + ^ng, no € Z, |no| < and k is the linear momentum k = ^n, n G 
[-L/2,L/2- 1] nZ. 

In the functional representation, the operator V of Eq.(2.1) is substituted with the following function 
of the Grassmann variables (3.2): 


V{ip,ip) = -h u>j j dt Y 




(3.6) 


je{i.2,3} 


—Ao 




(Q.a')G{t4}^ 

r 

/,+ I'n 


je{i, 2 , 3 } 


(0, (0, it)). 


Notice that V only depends on the fields located at the site a; = 0. This is important because it will 
allow us to reduce the problem to a 1-dimensional one [4,5]. 

The average of a physical observable F localized at a; = 0, which is a polynomial in the fields 
^^(0,t) and Pait)j will be denoted by 

{F)/=PJ P{dp)Po{d4,)e-^^^^^^ F, (3.7) 

in which Po{dip) is the Gaussian Grassmannian measure over the fields 3p^{0,t) localized at the site 0 
and with propagator 5 ^,q( 0 ,t) and Z is a normalization factor. 


2 


This means that all integrals will be defined and evaluated via the “Wick rule”. 



6 


The propagators can be split into scales by introducing a smooth cutoff function y which is different 
from 0 only on (j,l) and, denoting Np^=\og2l3, is such that = 1 for all \z\ € 

[f,Ni3]. Let 


m] \ dej 

gh, (0,^) = 2^ 


dkpdk ) 

(27r)2 —ikp + e(fc) 


x(2-^™((fc-cc7r/2)^ + fco")) 


(o>' t) - ffb (o> sir' (0) 


2=-N„ 


dr'(r) 'i' / f 


m— — N^ 


(3.8) 


where mo is an integer of order one (see below). 

Remark: The a; = ± label refers to the “quasi particle” momentum topp, where Pf is the Fermi 
momentum. The usual approach [8,17] is to decompose the field ij: into quasi-particle fields: 


V'a(0,t) = E 


(3.9) 


U) — -^ 


indeed, the introduction of quasi particles [8,17], is key to separating the oscillations on the Fermi scale 
from the propagators thus allowing a “naive” renormalization group analysis of fermionic models 
in which multiscale phenomena are important (as in the theory of the ground state of interacting 
fermions [8,9], or as in the Kondo model). In this case, however, since the fields are evaluated at a: = 0, 
such oscillations play no role, so we will not decompose the field. 

We set mo to be small enough (ie. negative enough) so that 2'^°pp < 1 and introduce a first 
approximation: we neglect and gE: replace e{k) in Eq.(3.5) by its first order Taylor expansion 

around ujpp, that is by ujk. As long as mp is small enough, for all m < mp the supports of the two 
functions x(2“^'"((A:—a;7r/2)^-|-fcg)), a; = ±1, which appear in the first of Eqs.(3.8) do not intersect, and 
approximating e{k) by uik is reasonable. We shall hereafter fix mp = 0 thus avoiding the introduction 
of a further length scale and keeping only two scales when no impurity is present. 

Since we are interested in the infrared properties of the system, we consider such approximations 
as minor and more of a simplification rather than an approximation, since the ultraviolet regime is 
expected to be trivial because of the discreteness of the model in the operator representation. 

After this approximation, the propagators of the model reduce to 




(27r)^ —ikp + ujk 


x(2-2™(fc2 + kl)) 



g-ifeo(T+0 ) 

-ikp 


xi^-^^kl). 


and satisfy the following scaling property: 

g5r'(0,T) = 2™5j'(0,2™r), gj^^r) = ffM(2-r). 
The Grassmannian fields are similarly decomposed into scales: 


0 0 

E 2t4-l±(0,2-“t), ptit)= t) 

m= — Njs m— — Njs 


(3.10) 


(3.11) 


(3.12) 
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with (0, t) and {t) being, respectively, assigned the following propagators: 

I = <5o.a'5i°l(2™(t - t'))- (3.13) 

Remark: by Eq.(3.11) this is equivalent to stating that the propagators associated with the (yjl™! 
fields are and respectively. 

Finally, we define 


4<™l±(0,t)"4^ ^ 2^4™'l±(0,t), E (3-14) 

m' — — N^ m' — — Nf 3 

Notice that the functions (C; t")) decay faster than any power as r tends to oo (as a conse¬ 

quence of the smoothness of the cut-off function y), so that at any fixed scale m < 0, fields (pl™l 
that are separated in time by more than 2“™ can be regarded as (almost) independent. 

The decomposition into scales allows us to express the quantities in Eq.(3.7) inductively (see (3.16)). 
For instance the partition function Z is given by 

0 

Z = exp (^ — P ^ (3.15) 

m——N0 

where, for Njs < to < 0, 

y[o](^[<o]^^[<o]) dg ^(^[<o]^^[<o]) 

in which S R and has no constant term, i.e. no fields independent term. 


4 Hierarchical Kondo model 

In this section, we define a hierarchical Kondo model, localized at x = 0 (the location of the impurity), 
inspired by the discussion in the previous section and the remark that the problem of the Kondo effect 
is reduced there to the evaluation of a functional integral over the fields 'ip{x,t),(p{t) with x = 0. The 
hierarchical model is a model that is represented using a functional integral, that shares a few features 
with the functional integral described in Sec.3, which are essential to the Kondo effect. Therein the 
fields and (z?!™! evaluated at x = 0 are assumed to be constant in t on scale 2“’”, to = 0, —1, —2 ..., 
and the propagators (0, r) and gjg^ (r) with large Matsubara momentum ko are neglected (< 7 ^“’'] = 0 
in Eq.(3.8)). 

The hierarchical Kondo model is defined by introducing a family of hierarchical fields and specifying 
a propagator for each pair of fields. The average of any monomial of fields is then computed using the 
Wick rule. 

As a preliminary step, we pave the time axis R with boxes of size 2“™ for every to € {0, —1,..., —Ng. 
To that end, we define the set of boxes on scale m as 

{[*2M,(* + l)2l-l)} (4.1) 

Given a box A G Q™, we define Ia as the center of A; conversely, given a point t G R, we define 
Z\[’"l(t) as the (unique) box on scale to that contains t. 



A naive approach would then be to define the hierarchical model in terms of the fields snd 

and neglect the propagators between fields in different boxes, but, as we will see below, such a 
model would be trivial (all propagators would vanish because of Fermi statistics). 

Instead, we further decompose each box into two half boxes: given A G Qm and rj G {—,+}, we 
define 

+ r;2-'"-2) (4.2) 

for m < 0 and similarly for to = 0. Thus Z\_ is the lower half of A and Z\+ the upper half. 

The elementary fields used to define the hierarchical Kondo model will be constant on each half-box 
and will be denoted by V'a"^^(^r;) and iy9[r^^(Z\^) for m G {0, —I,---, —Njs}, A G Qm, V & {—>+}) 
a e {t4}- 

We now define the propagators associated with ip and (f. The idea is to define propagators that 
are similar [18,19,11], in a sense made more precise below, to the non-hierarchical propagators defined 
in Eq.(3.4). Bearing that in mind, we compute the value of the non-hierarchical propagators between 

fields at the centers of two half boxes: given a box A G Qo and rj G let 6‘^=2~^ denote the 

distance between the centers of Z\_ and Z\+, we get 




E 


dkdkokosin{ko6) . , , 2 ,,def 
-F Wfc +ko)=va 


ivS) = V 


(27r 

dko sin(A:o^) 

27 r fcn 


*0 


def 


x(fco) = rib 


(4.3) 


in which a and b are constants, see [5, p.4465]. We define the hierarchical propagators, drawing 
inspiration from Eq.(4.3). In an effort to make computations more explicit, we set a = b = 1 and define 


((^M-(A_,)^M+(A,)) (4.4) 

for TO € {0, —1, • • •, —Np}, r] G {—, -l-}, A G Qm, a € {),, f}- All other propagators are 0. Note that if 
we had not defined the model using half boxes, all the propagators in Eq.(4.3) would vanish, and the 
model would be trivial. 


In order to link back to the non-hierarchical model, we define the following quantities: for all t G R, 

tf= E 2¥V'L”^l±(A[™+il(t)), E (4.5) 

m— — Nf3 m— — N^ 

(recall that to < 0 and Z\l™I(t) D (f)). The hierarchical model for the on-site Kondo effect 

so defined is such that the propagator on scale to between two fields vanishes unless both fields 
belong to the same box and, at the same time, to two different halves within that box. In addition, 
given t and t' that are such that \t — t'\ > 2“^, there exists one and only one scale that is 

such that Z\['"(*-‘')l(t) = and (t) ^ Therefore V(t,F) G R^, 

V(a,a') G {t4}^ 

{ip~{Q,t)ip+{0,t')) =(5Q,,a/2™<‘-‘')sign(t-t'). (4.6) 

The non-hierarchical analog of Eq.(4.6) is (we recall that ( • )^ was defined in Eq.(3.7)) 

0 

(V’a(0,t)e(0,O)K = <5«.«' E 2™'5fl(0,2-'(t-O) (4.7) 

from which we see that the hierarchical model boils down to neglecting the to' that are “wrong”, that 
is those that are different from mit-t'), and approximating by sign(t — t'). Similar considerations 
hold for ip. 



9 


The physical observables F considered here will be polynomials in the hierarchical fields; their 
averages, by analogy with Eq.(3.7), will be 


Zj 

(in which ( • ) is computed using the Wick rule and Eq.(4.4)) and, similarly to Eq.(3.6), 


(4.8) 


V{ip,‘p) = -h ^ iVj I dt (4.9) 

je{i.2,3} 


— An 


(a.a')e{t4}^ 


/ 1 O Q \ ^ 


je{i,2.3} 


in which tp^{0,t) and (p^it) are now defined in Eq.(4.5). 


Note that since the model defined above only involves field localized at the impurity site, that 
is at a; = 0, we only have to deal with 1-dimensional fermionic fields. This does not mean that the 
lattice supporting the electrons plays no role: on the contrary it will show up, and in an essential 
way, because the “dimension” of the electron field will be different from that of the impurity, as made 
already manifest by the factor 2™ 0 in Eq.(4.6). 


Clearly several properties of the non-hierarchical propagators, Eq.(3.10), are not reflected in Eq.(4.6). 
However it will be seen that even so simplified the model exhibits a “Hondo effect” in the sense outlined 
in Sec.l. 


5 Beta function for the partition function. 

In this section, we show how to compute the partition function Z of the hierarchical Hondo model 
(see Eq.(4.8)), and introduce the concept of a renormalization group flow in this context. We will first 
restrict the discussion to the h = 0 case, in which V = Vq; the case h ^ 0 is discussed in Sec.6. 

The computation is carried out in an inductive fashion by splitting the average (e'A'lV’.v) ^ into 
partial averages over the fields on scale m. Given m € {0, —1, • • •, — we define ( • )m the 

partial average over and for a G 4\ G Qm and rj G {—,+}, as well as 

and for Z\ G Q-m, rn < —Np, 

4^^”^1(Z\,)=^, (5.2) 

Notice that the fields and (p[r™~^^^(4\) play (temporarily) the role of external fields 

as they do not depend on the index ry, and are therefore independent of the half box in which the 
internal fields (4\^) and (/3[r™^^(4\^) are defined. In addition, by iterating Eq.(5.1), we can 

rewrite Eq.(4.5) as 

= = (5.3) 

We then define, for m G {0, —1, • • •, —Np}, 

y[o](^[<o]^^[<o]) i/o(^[<o],^[<o]) (5 4) 

in which G i? is a constant and contains no constant term. By a straightforward 

induction, we then find that Z is given again by Eq.(3.15) with the present definition of cl"*! (see 
Eq.(5.4)). 
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We will now prove by induction that the hierarchical Kondo model defined above is exactly solvable, 
in the sense that Eq.(5.4) can be written out explicitly as a finite system of equations. To that end it 
will be shown that can be parameterized by only four real numbers, a I™! = ■ ■ ■, 

and, in the process, the equation relating and (called the beta function) will be computed: 

3 

AGQm'iT'—^ r/—± 


where 

o'E(^)"= (5.6) 

in which and are vectors of polynomials in the fields, whose j-th component for j G 

{1,2,3} is 


aK™1^'(A) E 


(a,a')e{t.t}2 


5[<™]i(^) dg 


(a.a')e{T.t}2 



(5.7) 


For m = 0, by injecting Eq.(5.3) into Eq.(4.9), we find that can be written as in Eq.(5.5) with 
q:[° 1 = (Ao, 0,0,0). As follows from Eq.(5.13) below, for all initial conditions, the running couplings 
Q,[™] remain bounded, and are attracted by a sphere whose radius is independent of the initial data. 

We then compute using Eq.(5.4) and show it can be written as in Eq.(5.5). We first notice 

that the propagator in Eq.(4.4) is diagonal in A, and does not depend on the value of A, therefore, we 
can split the averaging over iIR\A±) for different A, as well as that over (pi™!(A). We thereby find 
that 

zi 

In addition, we rewrite 


rj^± 

rj—± k—0 n—0 

in which the sum over k only goes up to 2 as a consequence of the anticommutation relations, see 
lemma D.l; this also allows us to rewrite 

2 3 k ^ 

E M ^ (^)) = 1 + E (^) (5-10) 

k—O n—0 n—0 


where 


— CTn 


= a 


^2 — ^2 


= < 


— (^[("1)2 — 


(5.11) 
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At this point, we insert Eq.(5.10) into Eq.(5.9) and compute the average, which is a somewhat long 
computation, although finite (see App.B for the main shortcuts). We find that 


r}—-^ n—0 n=0 

with (in order to reduce the size of the following equation, we dropped all I™] from the right side) 

elm] ^ + 9(1 + 9el + 324£f 


c-Q T ~r ^<-2 \ tJ-i"±-C3 

^ (^0 + + 34^2 + 344 - 2£l 

^ 4 + 944 + T 


1 


1 CM V2' 


(5.13) 


4'"”'’= ^(24+ 3644+^^ 

1/1. 1. . 1^. 

24^° 


p[m-l\ 

''3 


Cl” 


2^:3 + 


The could then be reconstructed from Eq.(5.13) by inverting the map a 1—1 £ (see Eq.(5.11)). 

It is nevertheless convenient to work with the £’s as running couplings rather than with the ct’s. 

This concludes the proof of Eq.(5.5), and provides an explicit map, defined in Eq.(5.13) and which 
we denote by 7^, that is such that Einally, the appearing in Eq.(5.4) is given by 


p,lm] ^ _2^/3+'"log(C'[’”l) (5.14) 

which is well defined: it follows from Eq.(5.13) that > 1. 


The dynamical system defined by the map TZ in Eq.(5.13) admits a few non trivial fixed points. A 
numerical analysis shows that if the initial data Ag = oo is small and < 0 the flow converges to a fixed 
point £* 


4 = -a:o 


1 + 5x0 
1 - 4xo ’ 




where xg ~ 0.15878626704216... is the real root of 4 — 19x — 22x^ 
is (see Eq.(5.11)) 



- 107x3 = 


(5.15) 


0. The corresponding cx* 


= e, 


0> 


U) 


= £: 


2 ) 


— P* __ p*p* — _ 


12 


12 


^*2 

'0 


(5.16) 


Remark: Proving that the flow converges to £* analytically is complicated by the somewhat contrived 
expression of £*. It is however not difficult to prove that if the flow converges, then it must go to £* 
(see App.E). Since the numerical iterations of the flow converge quite clearly, we will not attempt a 
full proof of the convergence to the fixed point. 

Remark: A simpler case that can be treated analytically is that in which the irrelevant terms {£i and 
£ 3 ) are neglected (the flow in this case is (at least numerically) close to that of the full beta function 
in Eq.(5.13) projected onto = 4 = cst). Indeed the map reduces to 

Cl™! = l + Ul + 9£l 

= + - 2«?) (5-17) 


which can be shown to have 4 fixed points: 

/o = (0,0), unstable in the 4 direction and marginal in the 4 direction (repelling if 4 < 0,4 = 0), 
this is the trivial fixed point] 
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/+ = (0, |), stable in the £2 direction and marginal in the £q direction (repelling ii < 0,^2 = I), 
which we call the ferromagnetic fixed-point (because the flow converges to /+ in the ferromagnetic case, 
see below); 

/_ = (0, — ^) stable in both directions; 

f* = (—|, |), stable in both directions, which we call the anti-ferromagnetic fixed point (because the 
flow converges to f* in the anti-ferromagnetic case, see below). 

One can see by straightforward computations that the flow starting from — | < < 0 and ^ 2 °^ = 0 

converges to f* and that the flow starting from > 0 and ^ 2 ^ = 0 converges to /+ (see App.E). 


6 Beta function for the Kondo effect 

In this section, we discuss the Kondo effect in the hierarchical model: i.e. the phenomenon that as 
soon as the interaction is strictly repulsive {i.e. Aq < 0) the susceptibility of the impurity at zero 
temperature remains positive and finite, although it can become very large for small coupling. The 
problem will be rigorously reduced to the study of a dynamical system, extending the map i —>■ TZt in 
Eq.(5.13). The value of the susceptibility follows from the iterates of the map, as explained below. The 
computation will be performed numerically; a rigorous computer assisted analysis of the flow appears 
possible, but we have not attempted it because the results are very stable and clear. 

We introduce a magnetic field of amplitude h G R and direction u> G S 2 (in which ^2 denotes the 
2-sphere) acting on the impurity. As a consequence, the potential V becomes 

V{f;,(p) = Vo{i’,(p) - h f dt{ipt{t)(TYip~,{t))ujj (6.1) 

tefi.2.3}-^ 

The corresponding partition function is denoted by Zh^= {e ~^’') and the free energy of the system by 
fh^— — fi~^ log Zh. The impurity susceptibility is then defined as 

, .defd'^fh 

= (6.2) 

The h-dependent potential and the constant term, i.e. — and then defined in the same 

way as in Eq.(5.4), in terms of which. 


A= E 4r'- 

m— — N^ 


(6.3) 


We compute in the same way as in Sec.5. Because of the extra term in the potential in Eq.(6.1), 
the number of running coupling constants increases to nine: indeed we prove by induction that 
is parametrized by nine real numbers, = (q:q™^, • • •, Qfg")]) G i?®: 

^ ^ oK-l(Zi) (6.4) 

/iee™"=o »;£{-.+} 

where for n € {0,1, 2, 3} was defined in Eq.(5.6) and 

^K™](^)<f|/' lA[<n^](^).^ 

^K-](^)t/ iB[<-l(Z\).a; 

olGjn] i (A[^™1 (Z\) . u;) (b[^™ 1 (Z\) • u;) 


(6.5) 
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o'f^'CA) t/ i(B[^-l(A) . B[^™1(A)) (a[^- 1(A) • u;). 

We proceed as in Sec.5. For m = 0, we write Vh{ip,(p) as in Eq.(6.4) with cxh = (1,0,0,0, 0,/i, 0,0, 0). 
For m < 0, we rewrite 

(“p EE“irM-r'(/!)>,„=n( n Eh(E“So&"'(/s))\ m) 

A n,7] A r/—± k—0 n—0 

and, using lemma D.l, we rewrite 

E i = 1 + E41oLE’(^) (6-7) 

k—O n—0 n—0 


where is related to by Eq.(C.2). Inserting Eq.(6.7) into Eq.(6.6) the average is evaluated, 
although the computation is even longer than that in Sec.5, but can be performed easily using a 
computer (see App.I). The result of the computation is a map TZ which maps to as well 

as the expression for the constant c|("^. Their explicit expression is somewhat long, and is deferred to 
Eq.(C.l). 

By Eq.(5.14), we rewrite Eq.(6.2) as 


x{h,P) 







( 6 . 8 ) 


In addition, the derivatives of c|("^ can be computed exactly using the flow in Eq.(C.l): indeed dhCj^^ = 
and similarly for and can be computed inductively by deriving TZ{1.): 

) • 44“’, (6.9) 


and similarly for Therefore, using Eq.(C.I) and its derivatives, we can inductively compute 

X{l3,h). 

By a numerical study which produces results that are stable and clear we find that: 

(1) if Aq = oo < 0, aj = 0, j > 0, and h = 0, then the flow tends to a nontrivial, A-independent, fixed 
point £* (see Fig.6.1). 

We define nj{Xo) for j = 0,1,2,3 as the step of the flow at which the right-discrete derivative of 
ij/£* with respect to the step is largest. The reason for this definition is that, as Aq tends to 0, 
the flow of £j tends to a step function, so that for each component j the scale rij is a good measure of 
the number of iterations needed for that component to reach its fixed value. The Kondo temperature 
Pk is defined as and is the temperature at which the non-trivial fixed point is reached by all 

components. For small Aq, we find that (see Fig.F.I), for j = 0,1, 3, 

n,(Ao) =co|Aor^ + 0(l), cq « 0.5 (6.10) 

and (see Fig.F.2) 

n 2 (Ao) = C2|log2 |Ao||-bO(l), C2 « 2. (6.11) 


(2) In addition to the previously mentioned fixed point £*, there are at least three extra fixed points, 
located at .^4=^(6) 0; 0;0) and £^‘^='^(0,0, ±1/3,0) (see Fig.6.2). 
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Np 


Fig. 6.1 plot of ^ as a function of the iteration step Np for Ao = ao = —0.01. The relevant coupling £2 (curve 

number 2 , in green, color online) reaches its fixed point first, after which the marginal coupling £0 (number 0 , 
blue) tends to its fixed value, closely followed by the irrelevant couplings £i and £3 (number 1 , both are drawn 
in red since they are almost equal). 
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Fig. 6.2 phase diagram of the flow projected on the (4,4) plane, with initial conditions chosen in the plane 
that contains all four fixed points: £* (which is linearly stable and represented by a yellow circle), £0 (which 
has one linearly unstable direction and one quadratically marginal and is represented by a green cross), £+ 
(which has one linearly stable direction and one quadratically marginal and is represented by a red star), and 
£ 1 . (which is linearly stable, and is represented by a yellow circle). 


When the running coupling constants are at £*, the susceptibility remains hnite as /3 —> cx) and 
positive, whereas when they are at it grows linearly with /3 (which is why was called “trivial” 
in the introduction). 

In addition, when Aq < 0 the flow escapes along the unstable direction towards the neighborhood 
of £\, which is reached after n2(Ao) steps, but since it is marginally unstable for Aq < 0, it flows away 
towards £* after nic(Ao) steps. The susceptibility is therefore finite for Aq < 0 (see Fig.6.3 (which may 
be compared to the exact solution [7, Fig.3])). 

If Ao > 0, then the flow approaches £\. from the Aq > 0 side, which is marginally stable, so the flow 
never leaves the vicinity of and the susceptibility diverges as /3 —>■ 00. 

(3) We now discuss the flow at > 0 and address the question of continuity of the susceptibility in h 
as h —> 0. If Aq < 0 and ap = h 4 through £3 first behave similarly to the h = 0 

case and tend to the same fixed point £* and stay there until £4 through £§ become large enough, after 
which the flow tends to a fixed point in which £2 = 1/3, £5 = 2 and ij = 0 for j ^ 2,5 (see Fig.6.4). 

Setting the initial conditions for the flow as aj = a* for j = 0,1, 2,3 and as = h, we define rj{h) 
for j = 0,1,3 and j = 4, 5, 6, 7,8 as the step of the flow at which the discrete derivative of £j/£j is 
respectively smallest (that is most negative) and largest. Thus rj(h) measures when the flow leaves £*. 
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X 



Fig. 6.3 plot of x(/3i 0) as a function of log 2 P for Ao = —0.28. 


^C) ^‘1 ^3 



Nfi 

I log2('“)l 


Fig. 6.4 plot of ^ as a function of the iteration step Njs for Ao = —0.01 and h = 2 . Here £o through £3 

are the components of the non-trivial fixed point and IX through £.% are the values reached by £4 through 
Is of largest absolute value. The flow behaves similarly to that at h = 0 until £4 through £g become large, at 
which point the couplings decay to 0 , except for £5 and £ 2 . 


We find that (see Fig.F.3) for small h, 


rj{h) = Cr\og2h ^+0(1), Cr « 2.6. 


( 6 . 12 ) 


Note that the previous picture only holds if rj{h) log2(/3iy), that is PkIt- ^ 1- 

The susceptibility at 0 < ft. ^ is continuous in ft as ft —>■ 0 (see Fig.6.5). This, combined with 
the discussion in point (2) above, implies that the hierarchical Kondo model exhibits a Kondo effect. 

(4) In [20, Fig.17, p.836], there is a plot of as a function of For the sake of comparison, we 

have reproduced it for the hierarchical Kondo model (see Fig.6.6). 

Similarly to [20], we find that ^ seems to be affine as it approaches the Kondo temperature, 
although it is hard to tell for sure because of the scarcity of data points (by its very construction, 
the hierarchical Kondo model only admits inverse temperatures that are powers of 2 so the portion 
of Fig.6.6 that appears to be affine actually only contains three data points). However, we have found 
that such a diagram depends on Aq: indeed, by sampling values of |Ao| down to 10“"^, has been 

found to tend to 0 faster than {log but slower than {log ■ In order to get a more precise 

estimate on this exponent, one would need to consider |Ao| that are smaller than 10“^, which would 
give rise to numerical values larger than and since the numbers used to perform the numerical 

computations are yi86-extended precision floating point numbers, such values are too large. 
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X 



h 

Fig. 6.5 plot of x(/3i for ^ ^ 10“® at Ao = —0.28 and /3 = (so that the largest value for ph is ~ 1). 


Pk 

x(/3.0) 



Fig. 6.6 plot of as a function of ^ for various values of Aq: Aq = —0.024 (blue), Ao = —0.02412 (green), 

Ao = —0.05 (red). In [20], Ao = —0.024 and —0.02412. Note that the abscissa of the data points are 2“” for 
n > 0, so that there are only 3 points in the range [0.25,1]. The lines are drawn for visual aid. 


7 Concluding remarks 

(a) The hierarchical Kondo model defined in Sec.4 is a well defined statistical mechanics model, for 
which the partition function and correlation functions are unambiguously defined and finite as long 
as /3 is finite. In addition, since the magnetic susceptibility of the impurity can be rewritten as a 
correlation function: 

rP 

X{P,0)= dt{{ip+{0)cTip {0)){ip+{t)crip (7.1) 

do 

x(/3, 0) is a thermodynamical quantity of the model. 

(b) The qualitative behavior of the renormalization group flow is unchanged if all but the relevant 
and marginal running coupling constants (z.e. six constants out of nine) of the beta functions of Sec.5,6 
are neglected {i.e. set to 0 at every step of the iteration). In particular, we still find a Kondo effect. 

(c) In the hierarchical model defined in Sec.4, quantities other than the magnetic susceptibility of 
the impurity can be computed, although all observables must only involve fields localized at a; = 0. 
For instance, the response to a magnetic field acting on all sites of the fermionic chain as well as the 
impurity cannot be investigated in this model, since the sites of the chain with x ^ 0 are not accounted 
for. 

(c.a) We have attempted to extend the definition of the hierarchical model to allow observables 
on the sites of the chain at a; 0 by paving the space-time plane with square boxes (instead of paving 
the time axis with intervals, see Sec.4), defining hierarchical fields for each quarter box and postulating 
a propagator between them by analogy with the non-hierarchical model. The magnetic susceptibility 
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of the impurity is defined as the response to a magnetic field acting on every site of the chain and on 
the impurity, to which the susceptibility of the non-interacting chain is subtracted. We have found, 
iterating the flow numerically, that for such a model there is no Kondo effect, that is the impurity 
susceptibility diverges as f3 when /3 —)• oo. 

(c.b) A second approach has yielded better results, although it is not completely satisfactory. The 
idea is to incorporate the effect of the magnetic field h acting on the fermionic chain into the propagator 
of the non-hierarchical model, after which the potential V only depends on the site at x = 0, so that 
the hierarchical model can be defined in the same way as in Sec.4 but with an h-dependent propagator. 
In this model, we have found that there is a Kondo effect. 


A Comparison with the original Kondo model 


If the partition function for the original Kondo model in presence of a magnetic field h acting only on the 
impurity site and at finite L is denoted by Z%{l3, Ao, h) and the partition function for the model Eq.(2.1) with 
the same field h is denoted by ZK{f, Ao, h), then 



Zk{I3, Ao, h) — Zj({j3, Ao, h) + Zj({f3, 0 , 0) 

(A.l) 

so that by defining 

def .. Zff{l3,0,0) 

^ Z%{p,Xo,h) 

(A.2) 

we get 

mK{{3,Xo,h) = imK(/3, Ao,/i). 



mK{{3,Xo,h) = KmK{l3, Xo,h) 

(A.3) 


XK{l3,Xo,h) = ^XK{l3,Xo,h) + ^ l3mK{l3,Xo,h) 



xk{P,Xo,h) = KXK{l3,Xo,h) - {n - l)pmK{l3,Xo,hf. 



In addition 1 < k < 2: indeed the first inequality is trivial and the second follows from the variational principle 
(see [16, theorem 7.4.1, p.l88]): 

log Zf {13, Xo,h) = max(s(p) - fj.{Ho + V)) 

> s(po) - go{Ho) + MV) = siM - MHo) = log zf,{13, 0, 0) (A.4) 

where s(/i) is the entropy of the state p, and in which we used 

MV) = Tr V)IZk{I3, 0, 0) = 0. (A.5) 

Therefore, for f3hf ^ 1 (which implies that if there is a Kondo effect then Pm\ ^ !)> the model Eq.(2.1) 
exhibits a Kondo effect if and only if the original Kondo model does, therefore, for the purposes of this paper, 
both models are equivalent. 


B Some identities. 


In this appendix, we state three relations used to compute the flow equation Eq.(5.13), which follow from a 
patient algebraic meditation: 

(Af Af > = (2 + 

(AfAfAf) = (B.l) 

(AfAfA«A^/> =45,y 


where the lower case a denote (Ai ) = (A 2 ) and _ {ip+)■ 
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C Complete beta function 

The beta function for the flow described in Sec .6 is 

4 "*“" = + 1844 + 344 + 344 - 2fg) 


4’"“'' = ^(2^1 + 9^24 + 2^ + + 244 + + g44 + -/ o ) 

4 "*“" = ^(24 + 3644 + 4 + 64 + \tl + ^4 + 244 + ^ 44 ) 

4“-' = + i'-'* + 5'» + 55'»'' + ^4 + jlfeh + 

4"* ~ (j\m] 644 + 44 + 10844 + 1844 + 3ti4 + 64^7 + 44) 

= ^(24 + 1244 + 244 + 21644 + 644 + 36£i4 + 1244 + 244) 
^[m 1 ] _ — _l_ _l_ 3£g£g -|_ 344 + 774^5 + 1844 + 3t24 + 3fi4 

G 2 

+244) 

4"*“'' = ^(5^7 + + ^^44 + ^44 + ^44 + ^44 + ^44 

+ j^44) 

4™ = C^rn] 0^®^® ^44) 

(j{^\ + (4 + 4)" + 9f? + 9fi + 324£l + + 18^7 + 18^1 


in which we dropped the exponent on the right side. By considering the linearized flow equation (around 
4 = 0 ), we find that 4,4)454 are marginal, 4,4 relevant and irrelevant. The consequent linear 

flow is very different from the full flow discussed in Sec.6. 

The vector I is related to ol and via the following map: 

1 2 1 2 
4 = Qo, tl = Ol + 4 = 02 + 

„ 1 2 1 1 1 1 1 2 

4 = 03 + —Qq + yxOoOe + -OiQ 2 + -0408 + -Q5O7 + —Qe 

12 18 2 6 6 36 

+ ^000405 + ^'^ 1*^5 + ^0204 + —O4O5Q6 + 2^®'^^®'^® (C'- 2 ) 

4 = 04, 4 = 05, 4 = 06 + -O4O5 

n 1 1 1 1 2 

4 = 07 + -00O4 + -oiOs + -04O6 + —O4O5 

n 111 1 2 

4 = og + -0005 + -0204 + -05O6 + —04O5. 


D The algebra of the operators Ori,±- 

Lemma D.l Given 77 G { — ,+}, m <0 and A G Qm, the span of the operators {oky™'(Z\)}„g{o.i, 2 , 3 } defined 
in Eq.(5.6) is an algebra, that is all linear combinations of products of ol^ffi^ A) ’s is itself a linear eombination 
of {A) ’s. 

The same result holds for the span of the operators (Z\)}„g{o,. ., 8 } defined in Eq.(6.5). 

Proof. The only non-trivial part of this proof is to show that the product of two Onffs is a linear combination 
of On^fj S- 

Due to the anti-commutation of Grassmann variables, any linear combination of and squares 

to 0. Therefore, a straightforward computation shows that V(i,j) G {1,2,3}^, 

BliBli = 24.i+Gt+t 
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where the labels and (Zl) are dropped to alleviate the notation. In particular, this implies that any product 
of three for i £ {1,2,3} vanishes (because the product of the right side of the first of Eq.(D.l) and any 
Grassmann field vanishes) and similarly for the product of three 

Using Eq.(D.l), we prove that spanjoL^^'(Zl)}„g{o i 2 3 } is an algebra. For all n £ {0,1, 2, 3}, p £ {1, 2, 3}, 

«e{i,2}, 


Op^O, O3O„ = 0, OiOo = 0, Oo 


-O 3 , O 1 O 2 = 2*^3 


(D.2) 


(here the [-"‘I, (Zl) and ^ are dropped). This concludes the proof of the first claim. 

Next we prove that span{oL7,7*k-^)}ne{0. ■■,8} is an algebra. In addition to Eq.(D.2), we have, for all 

p£ {0,---,8}, 


O0O4 = xOt, O0O5 = xOs, OoOe = TxOs, O0O7 — OoOg = 0 , O1O5 = xO?, 

6 6 18 2 

O 1 O 4 = OiOe = O 1 O 7 = OiOs = 0, O 2 O 4 = 2 ^®’ O 2 O 5 = O 2 O 6 = O 2 O 7 — O 2 O 8 = 0, 

OaOp = 0, O 4 = qOi, O4O5 = 2 *^ 6 ’ O 4 O 8 = qOs, O4O7 = 0, O 5 = g02, (D-3) 

O5O7 = xOs, O 5 O 8 = 0, 06 = Tc*^ 3 , OeOr = 0608 = 0, O 7 = 0| = O 7 O 8 = 0. 

6 18 


This concludes the proof of the lemma. 


E Fixed points at ft, = 0 

We first compute the fixed points of Eq.(5.13) for £2 > 0. It follows from Eq.(5.13) that if € is a fixed point, 
then £1 = 6 £ 3 , which implies 

(1 - 3 ^ 2 ) (£ 2(1 + 3£2) + 6£? + £ 0 ) = 0. (E.l) 

If £2 > 0, Eq.(E.l) implies that either £2 = £1 = £0 = 0 or £2 = |. In the latter case, either £0 = £1 = 0 or 
£0 7 ^ 0 and Eq.(5.13) becomes 

f 3£g + 2£o + 6£i(3£i-l) = 0 


£i(l + 18£?)+£g(3£i-i) = 0. 


In particular, £i(l — 12£i) > 0, so that 


£0 = ±2 


£i(l + 18£{) 

1 - 12£i 


which we inject into Eq.(E.2) to find that £0 < 0 and 


(E.2) 


(E.3) 


1 - ^(3£i) + y (3£i)" - y (3£i)" + 107(3£i)^ = 0. 
Finally, we notice that is a solution of Eq.(E.4), which implies that 

4 - 19(3£i) - 22(3£i)^ - 107(3£i)^ 

which has a unique real solution. Finally, we find that if £1 satisfies Eq.(E.5), then 


£i(l + 18£{) _ l + 15£i 

1 - 12£i ^1 - 12£i■ 

We have therefore shown that, if £2 > 0, then Eq.(5.13) has three fixed points: 


£S - (0,0,0,0), £; := (^0,0,y0j , 


:= -xo 


1 + 5a;o Xo 1 xo 
1 - 4x0 ’ T’ 3’ 18 


(E.4) 


(E.5) 


(E.6) 


(E.7) 


In addition, it follows from Eq.(5.13) and Eq.(5.11) that, if Ao < 0, then (recall that a)}* = Ao and = 0, 
i = 1,2,3) 

£{,"*’ <0, 0 < £lr> < y 0 < £^’ < 6£[r’ < y (E.8) 
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for all m < 0, which implies that the set {£ \ £o < 0, £2 > 0, £1 > 0, £3 > 0} is stable under the flow. In 
addition, if > “fi then < £\^\ so that the flow cannot converge to £0 or £*^. Therefore if the flow 

converges, then it converges to £*. 

We now study the reduced flow Eq.(5.17), and prove that starting from —2/3 < < 0, ^ 2 °' = 0, the 

flow converges to /*. It follows from Eq.(5.17) that < 0, > 0 for all m < 0, so that if Eq.(5.17) 

converges to a fixed point, then it must converge to /*. In addition, by a straightforward induction, one finds 
that > ^ 2 ™^ if ^ 2 "' < Furthermore, (2£2™' + (^ 0 ”')^) < which implies that Therefore 

^ 2 "*' converges as m —^ — 00 . In addition, if £\^^ > and > — | — so that 

converges as well as m —> — 00 . The flow therefore tends to /*. 

Finally, we prove that starting from > 0,^2°^ = 0, the flow converges to /+. Similarly to the anti¬ 
ferromagnetic case, ^ 2 ™^ > 0 for all m < 0, < | and In addition, by a simple induction, 

if Ao < 1, then > 0 and + | “ ^ 2 ”*' i® strictly decreasing and positive. In conclusion, fg"' and ^ 2 "*' 
converge to /+. 


F Asymptotic behavior of nj{\o) and rj{h) 

In this appendix, we show plots to support the claims on the asymptotic behavior of nj{Xo) (see Eq.(6.10), 
Fig.F.l and Eq.(6.11), Fig.F.2) and rj{h) (see Eq.(6.12), Fig.F.3). The plots below have error bars which are 
due to the fact that nj(Ao) and rj{h) are integers, so their value could be off by ±1. 


™o(Ao)|Ao|, ni(Ao)|Ao|, nsfAollAo 



logiolAol 


Fig. F.l plot of nj(Ao)|Ao| for j — 0 (blue, color online) and j = 1,3 (red) as a function of | logj^g |Ao||. This 
plot confirms Eq.(6.10). 


712 (Aq) 

I log2 l-^oll 



Fig. F.2 plot of n 2 (Ao)| log 2 |Ao|| ^ as a function of | logj^g |Ao||. This plot confirms Eq.(6.11). 
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Fig. F.3 plot of rj{h) \ log 2 (/i)|| as a function of | log 2 (/i)|. This plot confirms Eq.(6.12). 


G Kondo effect, XY-model, free fermions 


In [1], given u £ [1, • • •, L], the Hamiltonian Hh = Ho—ha^, with 




n^n + 1 


+ 




(G.l) 


has been considered with suitable boundary conditions (see App.H), under which Ho and ctq + 1 are unitarily 
equivalent to (— cosq) a^a~ and, respectively, to 'jr ^ in which are fermionic creation 

and annihilation operators and the sums run over q’s that are such that = —1. It has been shown, [1]®, 
that, by defining 


2 h ^ 


Fn(C) = l + 9E 


L ^ C + cos q 

Q 

2 h 


F{z) = lim Fl{z) = 1+ — 


I 

Jo 


dq 


TV Jo (2 + COS q) 


(G.2) 


the partition function is equal to Z£('l in which Z£ is the partition function at h = 0 and is extensive (i.e. of 
0(6=°"-**^)) and (see App.H, Eq.(H.12)) 


log Cl (^.^) 


-I3h + 


2^£log(l + ®''^) 


. Pl{z) . 


dz 


(G.3) 


where the contour C is a closed curve which contains the zeros of Fl{C) {e.g. , for L —>■ 00 , a curve around the 
real interval [—1, y/l + if /j, < 0 and [—\/l + 1] if h > 0) but not around those of l + (which are on 

the imaginary axis and away from 0 by at least ^). In addition, it follows from a straightforward computation 

that {F{z) — Vj/h is equal to the analytical continuation of 2{z‘^ — from (l,oo) to C \ [—1,1]. 

At fixed P < 00 the partition function (l{P, h) has a non extensive limit C(/9, h) as L ^ 00 ; the C(/3, h) and 
the susceptibility and magnetization values m(/3, h) and h), are given in the thermodynamic limit by 


so that x(/l) 0) 


logC(^, h) 
m{P, h) 
xiP, h) 


01 PI 2/1 

-Ph+T—.f -;^l 0 g(H-r 

SvTi Jq 1 + (^2 — 1) 2 


-1+^ 


dz 


m Jc I + eP^ (^2 _ 1) 5+2/1 
2 C 1 dz 


■KiTcl + eP^ ((^2 _ i)i+2h)2 


2 sinh(/3) 
(l+cosh(/3)) 


and, in the /3 —>■ 00 limit. 


m(oo, h) 


2 h 

Vl + 4/i2 ’ 


x(oo,/i) 


2 

(1 +Y/i2)3/2 


(G.4) 


(G.5) 


^ see [1], Eq.(3.18) which, after integration by parts is equivalent to what follows. Since the scope of [1] 
was somewhat different we give here a complete self-contained account of the derivation of Eq.(G.2) and the 
following ones, see App.H. 
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both of which are finite. Adding an impurity at 0, with spin operators tq, the Hamiltonian 

H\ = Ho — h{al + Tq )—ActoT'o (G.6) 

is obtained. Does it exhibit a Kondo effect? 

Since tq commutes with the cr„ and, hence, with Hq, the average magnetization and susceptibility, 
/i, A) and /i. A), responding to a field h acting only on the site 0, can be expressed in terms 

of the functions ^(/3, h) and its derivatives CiP, h) and C'iP, h). By using the fact that ({p, h) and C'iPi s-re 
even in h, while is odd, we get: 


x‘”*(/?,0) = /3-'aIlogTr ^ (e 

T = ±l ^ 


— j3Ho+/3X(T^T+/3h((T^+T) 


p-^dilogZ^"^(p,h,A) 


(G.7) 


= /3- 


^ C' + C/3T + (C + prOPT (C + PtO''^^ 

/ ^ Vint \ / ^ 


^int 


= xiP, |A|) + P{m{P, |A|) + 1)^ + oo 


Since 0) is even in A, it diverges for P ^ oo independently of the sign of A, while x(/3, 0) is finite. Hence, 

the model yields Pauli’s paramagnetism, without a Kondo effect. 


Remarks: (1) Finally an analysis essentially identical to the above can be performed to study the model in 
Eq.(2.1) without impurity (and with or without spin) to check that the magnetic susceptibility to a field h 
acting only at a single site is finite: the result is the same as that of the XY model above: the single site 
susceptibility is finite and, up to a factor 2, given by the same formula x(/3, 0) = 


(2) The latter result makes clear both the essential roles for the Kondo effect of the spin and of the noncom¬ 
mutativity of the impurity spin components. 


H Some details on App.G 

The definition of Hh has to be supplemented by a boundary condition to give a meaning to (Tl+i- If = 
(a^ ± in'll) 12 define A/<n as n'la~ = and JV = JV<l- Then set as boundary condition 


it rfe/ j- 
^i + 1 


r 1 \AI ± 

= - (-1) CTl 


(parity-antiperiodic b.c.) so that Hh becomes 


1 

Hh = -h{2cytoP - 1) - - -|- cr“(-l)^"o-++i) 

^ n=l 

(-1)'^^ (-<^r ) + crZ (-1)'^^ ) )• 

Introducing the Pauli-Jordan transformation 

ih / I \-^<n i i zb 

ttn = (-1) (^n , ^L + l = • 


In these variables 


1 

Hh = -h{2aUP 


^n+1 ^n^n + l) 


Assume L =even and let l'^= {q\q = n = 0,1,..., I’ — 1}; then 


Hh = - L 

q q,q' 


+ A-g‘(q-qP'' 


- 1 ) 


A± ^ ^ = -1, qel 

In diagonal form let be a suitable unitary matrix such that 


Hh = Y.^. 




if = 


(H.l) 


(H.2) 


(H.3) 

(H.4) 


(H.5) 


(H.6) 
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Then Xj must satisfy 


hence 


(-^cosgA+A- - ^ |0> 

^ O n n' ' n" 


(A, +cosq)[/,,e-‘‘^‘' = 


| 0 > 


(H.7) 


Vg e /, where we used the fact that A~A^\Q) = 5p_q|0). We consider the two cases \j 7 ^ —cosq for all g e / 
or Xj = — cos qo for some qo € I- 
In the first case: 


N{Xj) Xj + cos q ’ 


provided FL{Xj)‘^= 1 + ^ 


2 h 


Xj + cos q 


= 0, 


(H.8) 


where N{Xj) is set in such a way that U is unitary, or, in the second case. 


Xj = — cos qo, Uj 


~ _y2 ^ 1 ,- 


that ( 


'Ui„i' = 0 . 


(H.9) 


Since — cosq takes values and the equation Fl{X) — 0 has ^ solutions, the spectrum of Hh is completely 
determined and given by the 2^ eigenvalues 

A(n) = y^njAj, n = (m,...,nz,), = 0,1 (H.IO) 

3 


and the partition function is 


log (/3, h)=Y, log(l + + y] log(l + ) = ^ log log(l + ). (H.ll) 

q>0 J jCI 

On the other hand, since the function F[(z)/F l{z) has L/2 poles with residue +1 (those corresponding to the 
zeros of Fl{z)) and L/2 poles with residue —1 (those corresponding to the poles of Fl{z)), the contour integral 
in the r.h.s. of Eq.(G.3) is equal to 

y] log(l + ) - Y, log(l + = y] log(l + ) - ^ log Zl (/3) = log Zz (/?, h) - log Zl (/3). (H.12) 

3 q >0 3 


I meankondo: a computer program to compute flow equations 

The computation of the flow equation Eq.(C.l) is quite long, but elementary, which makes it ideally suited 
for a computer. We therefore attach a program, called mecuikondo and written by I.Jauslin, used to carry 
it out (the computation has been checked independently by the other authors). One interesting feature of 
meankondo is that it has been designed in a model-agnostic way, that is, unlike its name might indicate, 
it is not specific to the Kondo model and can be used to compute and manipulate flow equations for a 
wide variety of fermionic hierarchical models. It may therefore be useful to anyone studying such models, so 
we have thoroughly documented its features and released the source code under an Apache 2.0 license. See 
http://ian.jauslin.org/software/meankondo for details. 
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